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Abstract 

We investigate the renormalization group flows and fixed point 
structure of many coupled minimal models. The models are coupled 
two by two by energy-energy couplings. We take the general approach 
where the bare couplings are all taken to be independent. New fixed 
points are found for models {N > 3). At these fixed points, the 
coupling constants all have the same magnitude, but some are positive 
while others are negative. By analogy with spin lattices, these can be 
interpreted as non-frustrated configurations with a maximal number 
of antiferromagnetic links. The stability of the different fixed points is 
studied. We compute the critical exponents and spin-spin correlation 
functions between different models. Our classification is shown to be 
complete. 
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A better understanding of coupled conformal field theories is primordial 
for the study of various problems in condensed matter theory as in the so- 
called spin-ladders or in disordered spin systems. Moreover, it could help to 
learn something about three dimensional theories. A promising approach to 
perturbed conformal field theory is the theory of integrable perturbations |jl|, 
^ ^. Here, we choose the more modest approach of perturbative conformal 
field theory introduced in the pionnering work of Zamolodchikov |Q. Namely, 
we study the fixed point structure of the renormalization group equations for 
several coupled identical minimal models. The couplings are initially taken 
to be independent, so that we have one dynamical variable for each pair of 
models. The RG equations are computed up the third power in the coupling 
constant by using the Coulomb gas representation of minimal conformal field 
theories. These techniques have already been used in order to study the 
effects of a weak disorder on second order phase transitions In this 

picture, the minimal models are parametrized by e, which is related to their 
central charge deviation from the Ising model value |. 

We present results for coupled models, > 3. For all of these models, 
we find symmetric fixed points (all the couplings are equal in magnitude and 
sign), and new "mixed" fixed points, that is, points where all couplings have 
the same magnitude, but different signs. The number and arrangement of the 
links is shown to be equivalent to a graph coloring problem. The solution to 
this problem is presented, and the resulting configurations are shown to be, 
by analogy to spin lattices, totally non-frustrated with a maximal number of 
antiferromagnetic links. The stability of the new fixed points is then studied, 
and we find that these are tricritical points. A study of the RG flows and the 
physical nature of the fixed points enables us to dress an accurate picture 
of the different phases. We then consider some of the critical exponents and 
the spin-spin correlation functions between points from different models. We 
end by looking at the physical significance of the new fixed points and by 
some comments on possible numerical verification of these results. 

We consider coupled minimal models with energy-energy couplings: 




Hq being the energy of a single unperturbed minimal model. The partition 
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function is given by 



N 



JV 




(2) 



.1=1 

and therefore can be seen as the perturbation of a conformal field theory of 
decoupled minimal models by the interaction term 



We will study the renormalization of the coupling constants. Assuming all 
couplings to be independent is not only the most general approach, but also 
the only reasonable one. 

For example, if we had begun with a chain of critical models coupled only 
to their nearest neighbours, first loop computations would show that the 
chain develops next nearest neighbour couplings. Higher order computations 
eventually couple every model to all others, in such a way that there are no 
perturbative fixed points compatible with the chain geometry. 

Renormalization group equations 

Let us first begin with the case = 3 depicted in Figure 1. 
To compute the renormalization of the coupling constants, we use the Coulomb 
gas representation of minimal models 0. In this picture, the central charge 
of the model is characterized by the parameter = | + e 



The parametrization of a+ is chosen is such a way that, for e equal to zero, 
we recover the Ising model. For instance, the three-states Potts model has 

The renormalization of gij is computed by expanding the interaction term, 
and then summing up the diagrams that contribute to ei{x)ej{x). The cal- 
culation techniques were used to study disordered systems (see for example 
P, and were found to be in good agreement with numerical simulations 
(at least for the three-states Potts model). Summing contributions up to two 
loops, and rewriting bare quantities in terms of renormalized ones, we get 
the following equations for three coupled models: 








e^l2 + ^13^23 - 2^12 (^13 + ^23) + 0(9^) 



= r 



dr 
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Piaisir)) = = egi3 + 912923 -^913 {gu + als) +^(9'^) 

/3(to(r)) = r^^ = e923 + 912913 -^923(9 u + 9n) +^(9^^) (5) 

Notice that we redefined (— 3e e) and {4:iTgij — > 9ij) in order to lighten 
notations. We are looking for fixed points of the renormalization group, that 
is, set of couplings such that f3{9ij) = for all This condition can be 
rewritten in a more convenient form as 

912 (^-e + ^(^13 + ^23)) = 913923 

(^-e+ ^(^?2 + ^13)^ = ^12^13, (6) 

Solutions are then easily obtained. Besides the trivial 9*^ = solution (stars 
denoting fixed point values), there is the symmetric solution: 

9*12 = 9*13 = 9*23 = i-^ + ^') + 0{e'), (7) 

which corresponds to the fixed point of N symmetrically coupled minimal 
models p. To see it, we project our space of coupling constants on the 
direction 9ij = g, and the beta functions reduce to 

r^ = e9 + 9'-9' + 0(9') (8) 

which has only one non-trivial fixed point 9* = (— e + e^) + 0{e^). 

The system (|]) also has an "anisotropic" or mixed solution given by 

913 = -73 913 = 923 = +73 (9) 

where 

l3 = e-e' + Oie^). (10) 

These are the only perturbative solutions (with no finite part for e = 0) 
modulo the permutation of models. We will generalize this result to the case 
of coupled minimal models. These calculations generalize easily to three 
different minimal models (parametrized by ei, £263) with similar conclusions. 
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The involved integrals have been performed in . Evidently, other solutions 
for models can be generated from solutions for a lower number of models 
by direct sum. 

Let us first comment the particular case e = 0, which corresponds to 
three coupled critical Ising models. If we consider three "isotropically" cou- 
pled Ising models, this system is described by a 3-colored Gross-Neveu model 
whose infrared behavior depends of the sign of the initial coupling con- 
stant (asymptotically free theory or massive regime) [|10]. When introducing 
anisotropy, the system (H) has a non-trivial fixed point (up to permutations) 
solution = fi'23 = which can be interpreted as one Ising model plus one 
Ashkin- Teller model. This fixed point can be reached only if g^2 < 0? fl'23 < 0- 
The 0(3) symmetry of the 3-colored Gross-Neveu model is then broken to 
SU{2) X U{1). 

We now return to the more general case e 7^ and study the stability 
of the new fixed points. It can be determined by looking at infinitesimal 
variations in the vinicity of the fixed points. If we take 

gi2 = 1 + Sgi2 gi2 = 1 + Sgi2 gn = -l + Sgi3, (11) 

then from the RG equations (|]), we have {g = {gu, 923, dis)) 

6t=ASg (12) 

where (it is sufficient to retain only first order in e) 



A 



e — e e 
— e e e 
e e e 



(13) 



Since e > (we redefined (e — 3e) earlier), the matrix A has two positive 
and one negative eigenvalues. The mixed fixed points are thus tricritical. 
The stable direction is given by 

Sgstab = {912 + 923 - 913) ■ (14) 

Calculations are easily generalized for an higher number of coupled mod- 
els. One still finds both isotropic and anisotropic solutions. If one defines 
the matrix g whose elements are gij, then the beta functions for general 
are 
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The condition (3{gij) = again leads to the symmetric solution. 

and to anisotropic configurations, with g*j = ±Jn, whose exact arrangement 
we will discuss below. Evidently, other solutions for N models can be gener- 
ated from solutions for a lower number of models by direct sum (a solution 
for N coupled models can be constructed with solutions for M and N — M 
models (M < N)). The value of the coupling 7jv is given by 

We can thus write all coTiplings as g*j = Tij^N with Tij = ±1. At this 
point, we can make an interesting analogy with lattice spin models. If Tij < 0, 
we can assume the coupling gij to be ferromagnetic (F) else antiferromagnetic 
(AF). If all bonds are ferromagnetic, Tij — —1. The existence of "mixed" 
solutions, that is solution with both F and AF couplings, was exphcitely 
shown for the case = 3 in a previous section. The anisotropic solution has 
two AF bonds and one F bond. Note that the product of the three coupling 
constants is also negative. By analogy to spin lattice systems (if we replace 
for example minimal models by genuine Ising spins) , it corresponds to a non- 
frustrated solution. Note that for Ising spins, this solution is equivalent to 
the symetric one. Finding the general configuration of ferromagnetic and 
antiferromagnetic bonds for N models is not obvious. However, by looking 
carefully at the renormalization equations, one can convince himself that 
solutions with couplings of different signs are possible if and only if every 
triplet of points i,j, k is such that 

TijTjkTik = -1. (18) 

which correspond to non-frustrated configurations. This is surely the case for 
the purely ferromagnetic solutions. To find other configurations, the problem 
can be reformulated in a more convenient way: 

Consider a system of N points. Points are linked two by two in such a 
way that every triplet of points contains one or three links. What are the 
possible configurations? Are they unique (up to permutations of models)? 
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This graph theory problem is identical to our RG fixed point condition 
if the presence of a link is the analogue of a ferromagnetic bond and the 
absence of link the analogue of an antiferromagnetic one. Besides the fully 
connected solution, there is a unique solution (up to permutation of models) 
such that every triplet of points is connected by one or three bonds. The 
total number of bonds, N^, is given by 

^"-{p^+p, N^2p+l. (^^^ 

The solution is constructed in the following way. Label sites from 1 to N. 
Then put links between sites a and [3 if and only if they both have even or 
odd indices. A given triplet will contain either indices of the same parity or 
two of one parity and one of the other. In the first case, there will be three 
links and in the other one. This solution has N^^ bonds. Proving the unicity 
of this construction is straightforward by induction. First, let us assume that 
there exists a chain of ferromagnetic links going from sites 1 to A'^ — 1 . Then 
there must be a ferromagnetic link between the site N and either site 1 or 2 
(not both). Both are equivalent by a permutation of models 1 and 2. The 
configuration is easily shown to be determined completely and leads to the 
above construction. Now, if there are no such chain between models 1 to 
A" — 1, we are in one of the two cases below: 

• There is a point with no ferromagnetic links. 

• There is a cluster of ferromagnetic links linking all models 

One can check that the second of these cases leads inevitably to the 
existence of a chain as in the above construction and thus to the same con- 
figurations. The first can only happen in a purely antiferromagnetic system. 
This establishes the unicity of the above construction. There are thus only 
two types of fixed points: 

• Isotropic ferromagnetic fixed point (each bond is ferromagnetic) 

• Anisotropic fixed points (arrangements of ferromagnetic and antifer- 
rromagnetic bonds) 

Having identified the different fixed points and their physical significance, 
we can now study the RG fiows and physical relevant quantities such as 
critical exponents and correlation functions. 
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To simplify calculations, we considered RG flows only for the case = 3. 
Higher dimensional behaviour is similar. We studied numerically the guQis 
projection of the flows for different values of (723- One can clearly view, for 
§23 = —73, both the ferromagnetic and mixed points. We have presented 
on Figure 2, the sketch of the flow for = 3 in the (5'i2,5'i3) plane. Fi 
and F2 are the two fixed point (plus the origin which is unstable in all the 
directions). A few remarks can be done. First, note that the "frustrated" 
configurations ((7i29'i35'23 > 0) correspond to the half plane gi2 > where 
the flow is always driven to strong couplings. Secondly, if we have the initial 
conditions gi2 < ^,913 > and \gi2\ > \gi3\, the flow is driven to strong 
coupling along the axis gi3 = 0. We recover a situation analogous to the case 
of three coupled Ising models. Namely, for e = | (the 3-states Potts model), 
we find in the infrared limit one critical Potts model plus two coupled Potts 
models which is known from integrable perturbation theory to be driven in 
a massive regime [0, 0. 

Critical exponents and correlation functions 

In this section, we will compute the correlation functions of energy and 
spin operators from different models. To do so, we must first compute the 
renormalized operators e' and a'. In the usual way, we are looking for ma- 
trices Ze and such that 



where a = (cti, cr2, ■ • • , (TAr),and similarly for e. A convenient way to define 
these constants is to consider an Hamiltonian of the form 



which can be seen as the original hamiltonian supplemented by source terms. 
The renormalization of the source terms is then done perturbatively by ex- 
panding the coupling term. 

Most of the involved calculations can be found in litterature |]^, the 
only difference here being that the source terms and couplings are assumed 
to be different for all models and pair of models. We will not repeat them 
here. Computing contributions and rewriting bare quantities in terms of 



(20) 
(21) 




(22) 
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renormalized ones, we obtain: 



dlog{Z,{r))ij _ 1,„ ^2 , ,n^„3 



gij-M' + Oig') (23) 



dr 2' 



dr 

with 



{1(9% (1 + ^) + ^{g'h + 0{g')^ 5,,, (24) 

^-V^(-1)P(-1)- ^''^ 

We are now able to compute the correlation functions. Prom the RG 
equations and conformal invariance of correlation functions, we have (for a 
lattice cut-off up to the scale ~ R) 

(£,(0)e,(i?)) ~ (26) 

and 

(a,(0)a,(i?)) r^Y,Y,{ZMZah-^MQ)ai{l)). (27) 

If we define 

= ^^M^ (28) 
dr 

and (7CT)ij similarly, we can replace the renormalized quantities by their fixed 
point value (denoted by a star). We then get 

(e,(0)£,(i?)) ~ J] J]i?-^^^°^-(-^)--(-^)^'(£,(0)eKl)) (29) 

and 

(a,(0)a,(i?)) - i?-^^^°^-(^^)--(^^)-(a,(0)a,(l)). (30) 

The simplifications for the spin-spin case are due to the fact that the renor- 
malization matrix Z„ is diagonal. For the energy-energy correlation func- 
tions, this is not the case and one must either diagonalize the matrix Z^ 
and introduce a new set of energy operators in order to define critical expo- 
nents or study long range behaviour where only the dominating term will be 
pertinent. 
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The shift of the critical exponent A^. is given by 

(A.),, = AW + i(7;, + 7a- (31) 

Computing 7*^ shows that the critical exponents for the anisotropic solutions 
are the same as the ones for the purely ferromagnetic case. This is the case 
because second order contributions contain the square of the couplings and 
that the third order terms obey the 

9ij9jk9ki = -1% (32) 

condition, the same for both fixed points. Any difference between the spin- 
spin critical exponent will be of order gf^ or higher. Things are different for 
the energy-energy correlation functions, and one finds different eigenvectors 
and eigenvalues in the purely ferromagnetic and mixed cases. The physical 
interpretation of the eigensystem is not obvious. 

As a conclusion, we give a brief summary of the results we found. By 
using perturbative conformal field theory techniques, we considered coupled 
minimal models in what we believe is the most general approach. This led 
us to the identification of new fixed points, which were shown, using graph 
theory arguments, to be, aside from the purely ferromagnetic solution, the 
only non-frustrated configurations. A study of the renormalization group 
fiows has shown that these new points are tricritical. Physically significant 
quantities such as critical exponents and correlation functions were com- 
puted perturbatively for the critical theories associated to the different fixed 
points. Spin-spin correlation functions were shown to be identical, up to g'^, 
for the ferromagnetic and mixed configurations. Energy-energy correlation 
functions however, turn out to be different and can be used to distinguish 
between both critical behaviours. It would be interesting to gain some nu- 
merical evidence of the existence of such fixed points for coupled 3-states 
Potts models. However, the tricritical nature of the mixed configTirations 
and the fact that couplings are asymetrical and therefore that the critical 
line is not self-dual certainly complicate computations. 
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CAPTIONS 

Figure 1: Three coupled minimal models represented by the big dots. 

Figure 2: A sketch of the projected renormalisation group flow for three 
coupled minimal models. 
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